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1. Einleitung 


Die Windenergie wird schon seit vielen hundert oder gar tausend Jahren 
genutzt. Die "Bockwindmühle" kannte man schon im Spätmittelalter. Das 
Mühlenhaus stand auf einem Bock (daher der Name) und wurde auf diesen g- 
dreht und der Windrichtung nachgeführt. Die Windmühlen konnten größer ge- 
baut werden, als man nur noch den Mühlenkopf der Windrichtung nachführte. 
Dieser Mühlentyp hei3t "Holländer Windmühle". Nun müssen Windräder aber 
nicht, wie diese, 4 Flügel haben. Die "Amerikanische Windturbine" aus dem 
19. Jahrhundert hatte oft mehr als 40 Flügel; der Monopterus, eine 
deutsche Forchungsanlage, hingegen hat nur einen Flügel. Die Windradachse 
kann auch vertikal angeordnet sein. Der Savonius- und der Darrieusrotor 
sind Beispiele für solche Vertikalachsrotoren. 


2. Der Wagner-Rotor 


In diese Reihe geschichtlicher Entwicklung und Verschiedenartigkeit glie- 
dert sich auch der Wagner-Rotor ein: Abb.i. Der Wagner-Rotor ist ein Win- 
drad mit einem Achswinkel 7 deutlich verschieden von O bzw. 90° und 
einem Konuswinkel € der Flügel deutlich verschieden von %0°. Der Wind 
strömt dabei in den, von den Flügeln gebildeten Kegel hinein. Da der Wag- 
ner-Rotor bei entsprechender Winkelwahl keinen Turm benötigt, und daher 
auch einen geringen Schwerpunkt hat, kann er auch auf einem Schiff mon- 
tiert werden und die hohen Windstärken auf See nutzen. 


Er Windradberechnung 


Bei der Windradberechnung geht man einerseits aus von den Grundgleichun- 
gen der Strömungsmechanik: Kontinuitätsgleichung, Energieerhaltungssatz, 
Impulssatz und andererseits von den gemessenen Größen Auftriebsbeiwert 
c_ und Gleitzahl € zur Beschreibung des Profilverhaltens unter Berück- 
sfchtigung vom Turbulenzgrad und der Zähigkeit. Die Strömung wird als in- 
kompressibel angesehen. Die Betrachtung ist zunächst eindimensional. Zu- 
sammen ergibt dies ein math./phys. Modell der Windradberechnung. Aus dem 
Modell erhält man sowohl Kenngrößen und Wirkungsgrade als auch Ausle- 
gungskriterien wie Flügeltiefe und Anstellwinkel. Die Berechnung wiri für 
einzelne Flügelelemente durchgeführt, die zusammen das Gesamtverhalten 
des Rotors ergeben. 
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4. Strömung am Flügelelement (siehe Abb. 2) 


Man kann sich leicht vorstellen, daß ein Flügelelement des Wagner-Rotors, 
das sich auf dem, von den Flügeln gebildeten, Trichter bzw. Kegelmantel 
bewegt, auf einer Kreisbahn um die Drehachse umläuft. Dieses Teilchen er- 
fährt durch die Drehung eine Geschwindigkeit tangential zur Kegelmantel- 
fläche, die Umfangsgeschwindigkeit u und durch die Windgeschwindigkeit 
eine Geschwindigkeit senkrecht zur Mantelfläche, genannt v. Man stellt 
sich nun das Flügelelement in Drehung an drei verschiedenen Orten vor: 

- vor dem Rotor in unbeeinflußter Strömung: Ebene | 

- in seiner wirklichen Stellung: Ebene 2 

- hinter dem Rotor, dort wo die Beeinflusssung der Strö- 

mung durch den Rotor abgeschlossen ist: Ebene 3 

Soll der Strömung Energie entzogen werden, muß die Geschwindigkeit v 
(Index steht für die Ebene) auf v, abgemindert werden. Wird das Windrad 
von der Strömung angetrieben, so wirkt nicht nur von der Luft ein Moment 
auf den Rotor, sondern auch ein Moment vom Rotor auf die Strömung. Dies 
verursacht einen Drall der Strömung und eine Vergrößerung der tangentia- 
len Geschwindigkeit von u, auf U,. Die Anströmgeschwindigkeit w ist die 
vektorielle Addition von u und v. 


5. Kenngrößen 


Aus dem math./phys. Modell leitet man zunächst Kenngrößen zur Beschrei- 
bung des Windrades ab: Die Schnellaufzahl A ist das Verhältnis aus Tan- 
gentialgeschwindigkeit u zur Geschwindigkeit senkrecht zum Kegelmantel 
V.» 
1 


A=suyı/vı Gl. 4 


Die effektive Schnellaufzahl A, beschreibt das gleiche Verhältnis, jedoch 
nicht in der Ebene 1, sondern in der Ebene 2, also am tatsächlichen Ort 
des Flügelelementes. 
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Die Blattspitzenschnellaufzahl A+ ist die Schnellaufzanl A , die an der 
Blattspitze auftritt. 

Der Abminderungsfaktor &£ beschreibt die Abbremsung der Windgeschwindig- 
keit bein Durchtritt durch den Rotor: 


E=v3/v Gl. 3 


Der Abminderungsfaktor kann Werte zwischen O und 1 annehmen. Eine Zahl 
nahe O0 bedeutet eine starke Abminderung. 

Der Drallfaktor 5 ist definiert als das Verhältnis der Tangentialge- 
schwindiskeiten: 


o=u3/uj Gl. 4 


Der Drallfaktor ist immer größer oder gleich 1. Weitere Kenngrößen sind 
der Auftriebsbeiwert , und die Gleitzahll & . 
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6. Wirkungsgradberechnung 


In der Windradberechnung kennt man drei Einzelwirkungsgrade und einen da- 
raus berechneten Gesamtwirkungsgrad. Aufbauend auf dem Modell zur Win- 
dradberechnung, trieben Leute wie Prandtl, Betz und Hütter die Windradbe- 
rechnung voran, so daß es in der Zeit um den 2. Weltkrieg gelang, Formeln 
für die theoretische Berechnung der Wirkungsgrade aufzustellen: 


6.1 Der ideale Leistungsbeiwert 
ist definiert als Quotient aus Windradleistung P und Gesamtleistung Po 
durch die Windradfläche. 


.=P 
BBi FO 3 
Mit "ideal" ist gemeint, daß ein Flügel ohne Verluste vorausgesetzt 
wird. Die theoretische Windradberechnung liefert die Formel 


=ı2 = 


o=Y1+(1-Z ] A Eis 7 


Der ideale Leistungsbeiwert ist also abhängig von der Schnellaufzahl und 
dem Abminderungsfaktor; Abb.3 zeigt dies. Bei unbeeinflußter Strömung, 
d.h.$=1, ist der leistungsbeiwert O0. Mit zunehmender Abminderung steigt 
er an, erreicht bei einer optimalen Abminderung 8 den größten Wert und 
fällt dann bei sehr starker Abminderung, wegen Es Iuftstaus und ver- 
mindertem Iuftdurchsatz, wieder ab. Weiterhin steigt der Leistungsbeiwert 
mit zunehmender Schnellaufzahl (siehe auch Abb.4). 

Der Leistungsbeiwert kann max. den Wert 16/27=0.593 annehmen. Dieser Wert 
wird bei einer Abminderung von 1/3 und unendlich großer Schnellaufzahl 
erreicht. Der Leistungsbeiwert Sn gilt zunächst nur für ein Flügelele- 
ment. 
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Abb. 3 Der ideale Leistungsbeiwert /1/ 
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6.2 Der Profilreibungs-Verlustfaktor 
ist folgendermaßen definiert: 


nützliche Vortriebsleistung am Flügel c1 8 


pP entnommene Leistung aus der Strömung 
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Er beschreibt also das Verhältnis Nutzen zu Aufwand am Flügel. Die theo- 
retische Windradberechnung liefert 


„Zile-ip) 
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Der Profilreibungs-Verlustfaktor ist also Abhängig von der 
Schnellaufzahl A , der Gleitzail &, der effektiven 
Schnellaufzahl A. und dem Abminderungsfaktor 8. Nimmt man eine optimale 
Abminderung an, so zeigt Abb.5 den Profilreibungs-Verlustfaktor aufgetra- 
gen über der Schnellaufzahl mit der Gleitzahl als Parameter. Sieht man 
einmal von sehr kleinen Schnellaufzahlen ab, so fällt der Profilreibungs- 
Verlustfaktor mit zunehmender Schnellaufzahl und abnehmender Gleitzahl. 
Auch der Profilreibungs-Verlustfaktor gilt zunächst nur für ein Flügele- 
lement. 


6.3 Der Blattzahlverlustfaktor 

berücksichtigt die induzierte Wirbelbildung an den Flügelenden. Diese 
Wirbel werden durch Ausgleichsvorgänge verursacht, die sich aus dem 
Druckunterschied am Flügel ergeben: Überdruck auf der Flügelunterseite 
und Unterdruck auf der Flügeloberseite. 

Der Blattzahlverlustfaktor berücksichtigt weiter die "luftstreifenbrei- 
te", also den Abstand der Flügel voneinander. Hat ein Rotor viele Flü- 
gel, so ist die Iuftstreifenbreite gering, die vom jeweiligen Flügel zu 
kontrollieren ist, und der Flügel kann seine Aufgabe besser erfüllen. 

Der Blattzahlverlustfaktor kann nicht exakt berechnet werden. Man ist 
hier auf Näherungsformeln angewiesen. Hütter, /1/, gibt eine solche For- 
nel an: 


n7=(1-1,39/[z/1#%,2))° sl. 10 


Diese Abhängigkeit des Blattzahlverlustfaktors von der Blattzahl und der 
Blattspitzenschnellaufzahl wird in Abb.6 wiedergegeben. Der Blattzahlver- 
lustfaktor gilt für den ganzen Rotor. 
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Abb. 4 idealer Leistungsbeiwert /2/ 
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Abb. 6 Blattzahlverlustiakior /1/ 
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6.4 Der Gesamtleis beiwert ist wieder definiert als Quotient 
aus Windradleistung P und Gesamtwindleistung Po durch die Windradfläche. 
Hierbei sind jetzt alle Verluste berücksichtigt. 


“pges =P/Py Gl. 19 


Die Gesamtwindleistung Po kann bei konstanter Windgeschwindigkeit nach 
der Formel 


Po=1/2°p*Az»v? 61, 12 


berechnet werden. Dabei ist P die Luftdichte und A, die Projektionsflä- 
che des Windrades zur Windrichtung. Ist der Gesamtfirkungsgrad bekannt, 
kann leicht die Gesamtleistung berechnet werden. Der Gesamtleisungsbei- 
wert gilt für den ganzen Rotor und berechnet sich, wie noch zu Zeigen 
ist, aus einer Integration bzw. Summation über der Rotorfläche mit den 
dort vorhandenen Einzelwirkungsgraden. 


6.5 Berechnung des Abminde faktors 

Bisher noch ungeklärt ist die Berechnung des Abminderunssfaktors g ‚der 
in die Berechnungsgleichungen der Einzelwirkungsgrade c. und 
Ye eingeht. Der Abminderungsfaktor kann nicht aus seiner DefinitP&n 


E=v3/viı | Gl. 3 


berechnet werden, da die Geschwindigkeit Vz nicht bekannt ist. Da die Ab- 
minderung letztlich ja vom Flügel abhängt, kann sie aus der Flügelgeome- 
trie berechnet werden: 
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Dabei ist B nur ein Platzhalter. Der Abminderungfaktor wird hier aus der 
Formel zur Flügeltiefe nach D. le Gourieres berechnet. /2/. Es muß fol- 
gendes beachtet werden: 
* Der Auftriebsbeiwert c_ ist vom Profilanströmwinkel d abhängig. 
* Der Profilanströmwink&l ist abhängig von der effektiven Schnel- 
laufzahl und dem Blattwinkelß 


a=arctan(1/AXg)-B Gl. 15 


* Die effektive Schnellaufzahl A, ist abhängis von der 
Schnellaufzahl A und eben wiederum von dem Abminderungsfaktor 8, 


_.‚1ito 
AorATze G1. 16 


mit© nach G1.7 . 
Eine Lösung für den Abminderungsfaktor ist also nur iterativ möglich. Da- 
bei setzt man zunächst A,=A und berechnet über d.und c_ eine erste Nähe- 
rung für& . Es Kann nun A,nach G1.16 und G1.7 genauör berechnet wer- 
den. 
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7. Programm zur Windradberechnung 


Das im Rahmen des Seminarvortrages angefertigte Program berechnet die 
über die Windradfläche gemittelten Wirkungsgrade, insbesondere den Ge- 
samtleistungsbeiwert Cpges - Das Programm berechnet alle Windräder, die 
der Beschreibung des Wagner-Rotors entsprechen (vgl. Abb. 1) und seine 
Grenzstellungen; damit z.B. auch den Horizontalachsrotor. 

Die Grobstruktur des Programms entnimmt man Abb.7 . Nach der Eingabe der 
Geometrie und den Berechnungsanweisungen wird in der ersten Schleife die 
unabhängige Variable, die Blattspitzenschnellaufzahl At, zur Berechnung 
der Kurven "Wirkungsgrade über Blattspitzenschnellaufzahl" vorgegeben. In 
den nächsten beiden Schleifen wird die Integration bzw. Summation über 
der Windradfläche, also über dem Umlaufwinkel © und dem Radius T, 
durchgeführt. Für jedes Flügelelement wiri zunächst, abhängig von seiner 
Lage (©,v) und der Blattspitzenschnellaufzahl die am Flügelelement 
wirksame Schnellaufzahl und der Abminderungsfaktor gemäß Abschnitt 6.5 
berechnet. Die anschließend berechneten Wirkungsgrade multipliziert mit 
der Fläche des Flügel-(bzw. Flächen)elementes werden jeweils summiert. 
Teilt man die Summen nach dem Durchlauf der © - Y-Schleifen durch die Ge- 
santfläche, so erhält man die gemittelten Wirkungsgrade "idealer Lei- 
stungbeiwert" und "Profilreibungs-Verlustfaktor". Der Blattzahlverlust- 
faktor 92 braucht nur einmal für jede Blattspitzenschnellaufzahl berech- 
net zu werden, da er für den gesamten Rotor gilt. Mit der Formel 


A 
“Pges = nz l/h 2 0p,°7p°44 61,17 


erhalten wir schließlich auch den Gesamtleistungsbeiwert Cpge- Ist die 
Berechnung für alle gewählten Blattspitzenschnellaufzanlen dürchgeführt, 
wird eine Wertetabelle und ein Plott ausgegeben. 


Programm zur 
Windradberechnu nr. 


| Schleile Ar 
i___ Schleife, S__) 
Berechnung A 
\terative Berechnun & 
Berechnung G, pi, 7p, AA, Ao 
| Summieren, jeweils 2: Co: aA 
7p- aA 
IPA 2 7p aA 
Berechnung 72 
A Ä 
Yepss= Aa Zip: paA Y 
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8. Ergebnisse 


Abb.3 zeigt die Kurven für einen Wagner-Rotor mit Achs- und Konuswinkel 
von je 45°. Man sieht, daß sich der Verlauf des Blattzahlverlustfaktors 
(durch einen Punkt "." gekennzeichnet) nicht von Abb.5 unterscheidet. Der 
ideale Leistungsbeiwert (+) steigt ähnlich Abb.3 an, fällt bei hohen 
Schnellaufzanlen aber wieder etwas ab. Für hohe Blattspitzenschnellauf- 
zahlen (über A« =6) ist dann der stark fallende Profilreibungs-Verlust- 
faktor begrenzend für den Gesamtleistungsbeiwert. 

Die Kurven für den Horizontalachsrotor (Abb.9) verlaufen ähnlich, der di- 
rekte Vergleich (Abb.10) zeigt jedoch, daß vor allem der Profilreibungs- 
Verlustfaktor und damit auch der Gesamtleistungsbeiwert beim Horizonta- 
lachsrotor höhere Werte liefert. 

Abb.11 zeigt die Auswirkungen der Vergrößerung der Achsneigung bei kon- 
stantem Konuswinkel. Der Gesamtleistungsbeiwert fällt mit zunehmendem 
Achswinkel. 

Abb.12 zeigt die Auswirkungen der Vergrößerung des Konuswinkels. Der Se- 
samtleistungsbeiwert steigt mit zunehmendem Konuswinkel. 

Die Ergebnisse des Programms haben gezeigt, daß die Wirkungsgrade des 
Wagner-Rotors sinken, je weiter man sich von den Winkeln des Horizonta- 
lachsrotors entfernt. Es muß weiter berücksichtigt werden, daß die Win- 
dradleistung P beim Wagner-Rotor nicht nur durch die schlechteren Wir- 
kungsgrade fällt, sondern auch dadurch, daß sich bei konstanter Flügel- 
länge die Projektionsfläche des Windrades zum Wind verkleinert. Trotzden 
hat der Wagner-Rotor seine Berechtigung bei der Nutzung der Windleistung 
auf See, wie eingangs dargestellt. An Land Können die Kosten für den Turm 
gegenüber dem Horizontalachsrotor reduziert werden. Ob dies die geringe- 
re Leistung bezogen auf die Flügsllänge ausgleicht, kann hier nicht ge- 
klärt werden. 


9. Zusammenfassung 


* P=P .c 
#6 REF aus der Mittelung von verschiedenen Wirkungsgraden 
üßet”ger Windradfläche berechnet. 
* Diese Wirkungsgrade erhält man aus einem math./phys. Modell, 
das 
Drall, 
induzierten Widerstand (Randwirbel), 
Profilverhalten, 
Reynoldszahleinflüsse (Turbulenzgrad), 
Reibung und Zähigkeit berücksichtigt. 
* Kompressibilität ist nicht berücksichtigt. 
* Der Wagner-Rotor schafft die Möglichkeit Windenergie 
auf See zu nutzen. 
* An Land ist der Wagner-Rotor dem Horizontalachsrotor 
unterlegen. - 
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10. Anhang 
10.1 Berechn der Schnellaufzahl A. am Flügelelement 





zZ I x 
| 





X 


Abb. 13 Koordinatenwahl an Wagner-Rotor 


Die Geschwindigkeit V_ ist die Windgeschwindigkeit. Es wird gefordert, 
daß der Wind nur von ifinen nach außen durch die Kegelfläche strömt. Diese 
Forderung ist erfüllt, wenn 
Gl. 19 
T<K 


gilt. Die Berechnung der Schnellaufzahl gliedert sich in drei Teile: 
Berechnung der Geschwindigkeit v. 
Berechnung der Geschwindigkeit u, .. 
Berechnung von A aus A=u, /v; i Me 
v.: Mit Hilfe der Differentialgeometrie wird die Komponente von v, sen- 
xhecht zur Kegelmantelfläche berechnet. Alle Angaben verstehen sch in 
Kart.-Koordinaten. Die Kegelmantelfläche wird beschrieben durch: 


Gl. 19 
x=(r»cos6, r*-sin®, r-eot«) OsxrsR 
D<Ox2e*1 
Die Windgeschwindigkeit wird beschrieben durch: 
et inz;0 Gl. 19 
Yo Vu SıNnT, F} ar . . a 


side 


Der Normalenvektor n steht senkrecht auf der Fläche des Kegelmantels. 


a ı — 
cos] |j-r»-sin® -cotk-cos® 

— [) 

n= }sinelX| r-cos8]| = r.|-cot«*»sino! Gl. 20a 
cot« 0 1 


. » [1 > . K3 » 
Die Projektion von v_ in Richtung ist 
j A 
-von 
w 


V 644.23 


wNn > 

In! 
Diese Gechwindigkeit senkrecht zum Kegelmantel ist die gesuchte Geschwin- 
digkeit v,. Es muß die Richtung -n gewählt werden, weil n nach innen in 
den Kegel gerichtet ist, von v_ aber die Komponente nach außen aus den 
Kegel hinaus gesucht ist. u 


cost-sint-cotk*co 
VjzV a=V „OST SINE COLE 70959 


WR De Gl. 21a 
Vvcot2«+1 


Auch die Umfangs- oder Tangentialgeschwindigkeit u, ist von der Windge- 
schwindigkeit abhängig. Sie vergrößert sich, wenn der Flügel dem Wind 
entgegen läuft und umgekehrt. 


y ur 


Be T-sno 
! ne 


we Boa Vw sn 





uf =6*r=wer Bi: 29 
U]=w«r-v «sint-sin® Gl. 23 


Es kann jetzt die Schnellaufzahl A aus den errechneten Geschwindigkei- 
ten ermittelt werden. 


‚.Ur,lwer-v .sint-sina)-Vcot?«+1 SET 
v v leost-sint-cot«-cos0) : 


BR Ep 


Die Blattspitzenschnellaufzahl am Wagner-Rotor wird zweckmäßig so defi- 
niert: 


1, -u1-urR Gl. 25 
1 Vi 


Da dimensionslose Kennzahlen berechnet werden sollen, ist die absolute 
Größe des Windrades nicht von Bedeutung. Darum wird der Maximalradius R 
einfach gleich 1 gesetzt: 
Re] Gl. 26 
Damit ist 
wW>= A .y Gl . 2 7 
t w 


Setzt man w nach G1.27 in G1.24 ein und kürzt die Größe V„ 50 erhält man 
die Indgleichung 

(A,er-sint-sine)-Voot?«+1 

= POSTER EIN OR PR REREEREANE RER SENEN. Gl. 23 


cost-sint-cot«-cos® 


Es muß beachtet werden, daß r nur Werte zwischen O und 1 annehmen darf. r 
gibt also das Verhältnis an aktudler Radius:Maximalradius beim 
untersuchten Windrad. 


an 





40 


10.2 Programmliste des Programms zur Windradberechnung mit Hinweisen für 
ogrammbenutzer 


20= PROGRAM CPw (INPUT,OUTPUTsTAPES=INPUT»TAPE6=OUTPUT) 
30=C 
LT EI ou 2. 272° 2 22:2 2 2.220202 202. 2.2.2722. 2.2°272 202.2. 27202.2 0202.22 027200020587 
s0=sC tus ur 
60=C ### PROGRAMM ZUR WINDRADBERECHNUNG #2% 
TO=SC #us wir 
8O=SC HunRuHa urn TE SE A 
90=C 
100= REAL LAMMINJLAMMAX,LAMT,LAMsLAMO,L,LMIN»KAP 
110= CHARACTER PL*1,OBPL#3,WASPL#4, ZAHLPL#=B,PARAM#=7,KTR#3 
120= DIMENSION CA(-20:31),EPSt-20:°91),PL(35s91)»ZAHLPL (35) „PAWERTI{100) 
130= COT(PHTI=L/TAN(PHI) / 
140= RAD(PHI)=PHI/LA0#PL | 
150= GRAD (PHI)=PHI*180/PI} 
160= PI=3.1415927 
170=C \ v4 
180=C *#» EINGABE: GEOMETRIE. B WEISUNGEN #»# 
190=C 
200= DATA LAMMINsLANMAX,DELLAM /1»205,1/ 
210= DATA KAP,TAU,DELTHE /90,45,20/ 
220= DATA L+LMIS /10»37 
230= DATA TTsTL4MIN /70.5,0.5/ 
240= DATA BETT,RETL“IN /0»0/ 
250= DATA Z /27/ 
260= DATA DELR„SENAU 70.050.017 
eT0= DATA KTR,NSPLs vASPL /FOFFr, IOFFr,11000:7/ 
280= DATA DPARAM /'Kap / 
e90= DATA PAXERT /90,67.5»45,97%*0/ 
300= DATA DAANZ 73/7 
310= DATA ca /E1E Os leo2se3s hy Dse6g. Tg, Bne9s1le09sl.1lr1.2,1.3H1l.2H1l.ls 
320= 1 1.9s.2,23#.3>15#.95102.8,57,.7,5%.6,5# ,555%.4,6%,2/ 
330= DATA EPS/21=0s10»,20»25%,30+33,40,45,45,40,30,25,20,8*#10, 
340= 1 19*2s19*1.3»10=0.9s10#*0.6»,10*9.4»10*0.2-+11*0.1/ 
350= DATA aLPMIiI 7/97 
360=C 
370=C *### 1. FINGSABEXONT-OLLE #rr 
3806=C 
390= IF(LAMMIN.CT.LAMMAX) CALL ERRORSC!LAMMIN.GT.LAMMAX I ,#5) 
400= IF(LAMMIN.LE.O) CALL ERROR6SCFLAMMIN.LE,0',#5) 
410= IF(DELLAM.LE.O) CALL ERRORSL(IDELLAM.LE.0',*5) 
420= IF(DELTHE.LE.O) CALL ERRORSCHDELTHE.LE.0',#5) 
430= IF(3EN/DELTHE. IE.INT(3607/DELTHE) J)DELTHE=360/( INTL360 /DELTHE) +]) 
440= IFtPAFTAM.NE. KAP „AND. PARAM.NE. !TAUt.AND.PARAM.NE, !L® 
450= 1 „AND.PArAM.NE,*'LMIN!® .AND-PARAM.NE. ıTTe® «AND.PARAM« NE. "TLMIN! 
460= 2 sANN.PARAM.NE.IBETTF JAND.PARAM.NE. T3BETTLMIN® JAND.PARAM.NEs Zr) 
470= 3. THes 
480= PAANZ=]L 
490= PAWERT(l)=N | 
500= END IF 
510= TAU=SRAN (Tat!) 
520= KAP=RAN(KAP) 


au 


5s30= DELTHE=RAND(DELTHE) 

S40=C_ 

550=C ### PLOTTVORBEREITUNG ##% 

560=C 

570= DATA ZAHLPL /'9.0 0.0'r,5* ı,r0.12 D.2r,6*8 tr 
580= 1 ın.24 0.4'.,.6#° t,10,.36 0.9! ,6%** tr 
590= 2 r),.47 0.8',6%#* r,10,59 1.08/ 
600= DATA PL /31Aa5#r 7 

610= DOSK=1,35 

620= a PL(K«1l)='!ft 

630= DO9IK=L1,91 

640= 9 PLtl»x)=!t-! 

650= IFtOBPL.NE. OFF*) DELLAM=(LAMMAX-LAMMIN)/9O 

660=C 

6T7T0O=C #e# AUSGARE GEHAENLTER GROESSEN ##* 

680=C 

690= WRITE (6»100) DELR»GENAU.GRAD (DELTHE) sKTR»PARAM» 
700= 1 (PAWERT(I) ,„I=1»PAANZ) 

710= 100 FORMAT (//7///1H »10XrFDELR : 1,F10.6/ 

T20= 1 IH »10Xs '!GENAU = °s,F10.6/ 

730= 2 l= »10X»+'DELTHE = 1,F10.6// 

T40= 3 IH »10Xr '!KTR 8 ,A3// 

750= 4 14 s10Xs "PARAMETER: 1,A7/ 

760= 5 14 »s10X» WERTE : 2. 10F10.3/ 

7T0= 6 914 »10X rt t,10F10.3/)) 

780=C 

790=C *### VERAENDERN DES PARAMETERS #*% 

800=C 

8l0= DO4II=1»PAANZ 

820= IF(PARAM.EN. KAP v) KAP=RAD(PAWERT(ITD)) 

830= IF(PARAM,.EQ,„ Tau ’) TAU=SRAD(PAWERT(II)) 

840= IF(PARAM,.EN. HL v) L=PAWERT(II) 

850= IF(PARAM,.EQ.'LHIN t) LMIN=PAWERT ID) 

860= IF(PARAM.EQ.FTT ?) TT=PAWERT (II) 

870= IFtPARAM.EN. TLMIN t) TLMIN=PAWERT (II) 

880= IF(PARAM.EQ. BETT ı) BETT=PAWERT (II) 

890= IFtPARAM.EN. YBETLMIN!) 3ETLMIN=PAWERT (II) 

900= IF(PARAHLFEN. Z t) Z=PAWERT (II) 

910=C 

920=C #au 2, EINGABEKONTROLLE ##* 

930=C 

940= IF(TAU.LT,.N CALL ERRORBEL'ITAU.LT.O',#5) 

950= IF(GRAN (KA) .5GT.30) CALL ERRORE(IKAP.GT.90r,#5) 
960= IF(Tay,ST,raPp) CALL ERRÜRE ( TauU,5T.KAPı,#5) 
970= IF(LMIN.LT,O) CALL ERRÜRE (HLMIN.LT.0',#5) 
980= IF(LMIN.ST.L) CALL ERROROLILMIN.ST.L',#5) 
990= IFeTTELTEN.ORTLUINELT.N).CALL ERROR6LITT BZwe TLMIN.LT.9',#5) 
1000= IF(Z.LT.0) Z=-Z 

1010= Z=INT(Z) 

1020= IF{TAUJER.I.AND.CRAD (KAP) „EQ.90) DELTHE=SRAD (350) 
1030= WRITE (&»2ON)GRaDIKAP) »GRADLTAU) sL»LMIN»TTSTLMINsSETTSBETLMINsZ 
1040= 2H90 FORMAT(IHl»10Xert Kap TAU L LMIN TT TLMIN . 
1050= 1 I3aeETT ZETLMIN  ze//lMn »10Xs9F7.3/////1H »12% 
1060= 2 tLAMT cPI ETaAP? ETA ETAZ CPGES'!/Y) 
1070=C 

1080=C *#= VORBEPFITUNG Hirt 

1090=C 

1100= RMIN=LMIN/L 

1110= Q=(BETLMIN-BETTIZEI-RMIN) 

1120= RADIUS=L*SIN(KAP) 

1130= T=TT/RADIUS 

1140= U= (TLMIN-TT)/C1I-RMIN)Z/RADIUS 

1150= J=0 

1160=C 

1170=C *## BERECHNUNG DER WIRKUNGSGRADE #*## 

1180=C 


a5 





1190= 
1200= 
1210= 


" 1220= 


1230= 
1240= 
1250= 
1260= 
1270= 
1280=C 
1290= 
1300= 
1310=C 
1320= 
1330= 
1340= 
1350= 
1360= 
1370= 
1380= 
1390= 
1400= 
1410= 
1420= 
1430= 
1440= 
1450= 
1460= 
1470=C 
1480= 
1490= 
1500= 
1510= 
1520= 
1530= 
1540= 
1550= 
1560= 
1570= 
1580= 
1590= 
1600= 1 
1610=C 
1620= 2 
1630=C 
1640= 
1650= 
1660= 
1670= 
1680= 
1690= 
1700= 
1710= 
1720= 
1730= 
1740= 
1750= 
1760= 
1770= 
1780= 
1790= 
1800= 
1810= 
1820= 
1830= 
1840=C 


DN3 LAMT=ZSLAMMIN„LAMMAX,DELLAM 
MARKE]=O 

MARKE?2=O 

MARKE3=O 

MARKE4=0 

J=y*+l 

CPISUM=O 

ETPSUM=Q 

FETASUM=O 


DO2ZTHE=0,2=PI-DELTHEsDELTHE 
THEM=THE+NELTHE/2 


DOl R=ERMINs1-DELR»DELR 

MARKEI1=0 

RM=ZR+DELR/Z 

LAM= (LAMT#RM=-SIN(TAU)#SINLCTHEM) ) #SQRT (COT (KAPI#=2+1)/ 
(COS(TAU)=-SINLTAU)#COT (KAP) #COS (THEM) ) 


CALL SUBXILLAM,CA»ALPMINsBETT,Q,T» Us ZARMSGENAUsÄI»ALP»IALP»LAMO, 


KTRsMARKEL»MARKEII»MARKE3,MARKE4,#3) 

SIG=SQRT (l+(1-xI##2) /LAM*F*2) 
CPI=LAM=#2*(]+XI)*(SIG-]) 

IF(CPI1I.65T.0.593) CALL ERRÜR2 (#3) 

IFtCPI.LT.O) CaLL ERROR4 (#3) 

IF(MARKEIL.EQ.1) CPI=-CPI 
EPSM=EPS (IALP)+(ALP-IALPI*(EPS(IALP+L)-EPSC(CIALP)) 
ETAP=2#LAM#F(EPSM-LAMO)/LI+ÄI)/ULL+EPSM#=LAMO) 
IF(ETAP.GT.1ITHEN 

## MARKED ## 

IF(KTR.ER. OFF) THEN 

ETAP=] 

MARKE2=? 

ELSE 

CALL ERROR3(#3) 

END IF 
END IF 

IFtETAP.LT.O) ETAP=O 
DELAK=0,5#SURT(COT(KAF)R#2+L)#((R+DELR) #H2-R##2)#DELTHE 
CPISUM=SCPISUM+CPI*#FDELAK 
ETPSUM=FTPSUM+ETAP=DELAK 
ETASUM=ETASUM+CPI*#ETAP*=DELAK 
CONTINUFS 


CONTINUE 


AKZSIRTICHT(XAP)F=2+ 1) F1I®PI 
ETA=SETASUMZAH 

IF(ETA.ST.1) CILL ERRORZ3(HB) 
IF(ETA.LT.)) TSull ERRORA (HZ) 
ETaZ=(1-1.39/(7#35GRTC1+LAMTFRZ2)))HR2 
IF(ETaZ.G6T.1) FALL ERRUR3 (#3) 
I=SNINTI34*ETAZ+]) ® 
IFtwWASPL(a:a) „EN. LI) PL(Isy)=t,t 
ETAP=ETPSUMZAK 

IF(ETaP.6ST.1) CALL ERROR3(*3) 
I=NINT(34#ETAP+]) 
IF(AASPL(3:3).ZQ. 1") PLtIsJI=t,t 
CPI=SCPISUM/AK 

IF(CP1.5T.0.593) CaLL ErROR2(#3) 
IF(CPI.LT.0O) CaLL ERRÜR4(#3) 
I=NINT(CPI/0.593*34+1) 
IF(wASPL (2:2) .E9. tl!) PL(I»J)=strt 
CPGES=ETA=ETAZ 
I=NINT(CPSES/.593*34+]) 
IF(WwASPL(1:1).ZQ. 1!) PLtIsJ)I=ta0 


a6 





1850=C *## AUSGAZE DER WIRKUNGSGRADE DES GESAMTEN ROTORS *### 
1860=C 


1870= WRITE (8»30°) LAMTsCPI»ETAP»ETA»ETAZsCPGES»MARKEI»MARKE2, 
1880= 1 MARKE3»MARKE4 

1890= 300 FORMAT(IMH s10Xs6F7.3»604,413) 

1900= IFtUEQ-SN.OR.J.EQ.110.0R.J»ER.170.0OR.J.ER.230) WRITE(A,301) 
1910= 301 FORMAT (IHN) 

1920=C 

1930= 3 CONTINUE 

1940=C 

1950= 4 CONTINUE 

1960=C 

1970=C *## PLOTTAUSGARE *r* 

1980=C 

1990= IFtORPL.EQ.FOFF*) GOTO 5 

2000= D= (LAMMAX-LAMMINIYS 

2010= WRITE(6»400) (ZAHLPL(DI »(PL(IsJ)sJ=1»91),1=35H 1, -1) >» 

2020= 1 LAMMINJLAMMIN+DsLAMMIN+2#D,LAMMIN+3#D,LAMMIN+4 RD ,LAMMAX 


2030= 490 FORMAT(IH1,1Xs #1, 117Xs r#t ////// 7/18 sl0X,! #r+ ..'/ 


2040= ı 35(1H »i0X»48s9lal/)/IM s16Xs6(F4.1»14%X)//// 
2050= e IH »„1Xs #0, 117%, 1#r 7///1H s10Xs ICPGES: #'/ 

2060= 3 IH »10X, "CPI +97 IH »10X, !ETAP 2 »r/ 

2070= 4 IH »10Xs HETAZ : .r7///) 

2080=C 

2090= s SToP 

2100= END 

2e1l0=C 

2120=C ### UNTERPROGRAMM: BERECHNUNG VON XI» ALP» LAMO ### 

2130=C 

2140= - SURROUTINE SUBXILLAMsCAsJALPMINJBETTSQSTsUrZERM»GENAU,XI» 
2150= ı ALPsIALPsLAMOsKTR,MARKEL»MARKELISMARKEI,MARKEA,#) 
2160= REAL LAM»LAMO 

2170= CHARACTER «TR#23 

2180= DIMENSION CA(-20:90) 

2190= GRAD LPHAI)=PHI*130/P1 

2200°= PI=3,.1415927 

2210= LAMO=LAM 

2220= AT=zU#(]-Rit)+T 

2230= xIA=0 

PZ40= G0TO2 

2250=C 

2260>= 1 Lamg=snRrT(##tLam##2+l)/(1+x1)#®2-1) 

22T0= xIA=AI 

2?280= 2 3ET=I8 (1) SM) #DETT 

2290= ALP=GPADLATANLLZ/ZLAMO) I) -RET 

2300= IF<ESWSLTSCH ALFA Pr+)20 

2310>= IFe ALP .LT.ALPMIN ORAL? .GT.90) .AND.KTR.NE.rOFFr) CALL ERROR! (#3) 
2320=C az MARKEL ## 
2330= TEeLSLPeLTRALDHINORLALP.GTEIO) .AND.KTR.EQ. OFF) THEN 
2340= IF(LaMO.LT.O) Lamn=-LAr0 

2?2350= IFtALP.LT.2#ALPMIN-90) „CALL ERRORL(3) 

2360= IF(aLP.LT.ALP: IV) ALPNZU=2#HALPMIN-ALP 

2379= IFtALP.GT.SÜ.AHn ,ALP.LT.130-ALPMIN) ALPNEU=]130I-ALP 
2380= IF(ALP.GT.130-aLFfH IN) CALL ERRORLI=3) 

2390= ALP=ALPNE! 

2400= MARKE]=S] 

24 10= MARKZI1=1 

2420= END IF 

2430=C 

2440= IALP=INT(ALP) 

2450= CAaa=sCat ALP) +LALP-IALP)=(CALIALP+L)-CALTALPI) 

2460= IFtCAA.LT.9I CALL ERRORS (#3) 

2+T0=C ## MARKEI3 ## 
?480= IFICAM.LT.LE-3) THEN 

2490= xI=1 

2500= MARKT3=3 


ar 





2510= 
2520= 
2530= 
2540= 
2550= 
2560=C 
2570= 
2580= 
2590= 
2600= 
2610= 
2620= 
2630= 
2640= 
2650= 
2660=C 
26T0O=C 
2680=C 
2690=C 
2700= 
2710= 
2T20= 
2730= 
2740= 
2750=C 
2760= 
2T7T0= 
2780= 
2790= 
2800= 
281l0=C 
2820= 
2830= 
2840= 
2850= 
e8B60= 
2ea70=C 
2880= 
2890= 
2900= 
2910= 
2920= 
2930=C 
2940= 
2950= 
2960= 
2970= 
2980= 
2990=C 
3000= 
3010= 
3020= 
3030= 
3040= 
3050=C 
3060=C 
3070=C 
3080=C 
3090=C 
3100=C 
3110=C 
3120=C 
3130=C 
3140=C 
3150=C 
3160=C 


RETURN 
EnD IF 
A=SBHPI/LCAMRZ) /QT#RM 
B=4/ (LAMO=SGRT (LAMOF#FZ+L)) 
XI=(8-1)/(3+1) 
## MARKE4 #*% 
IF(X1.LT.0) THEN 
XI=0 
MARKE4=4 
RETURN 
enD IF 
IF(ABSC(XIAa-XIY/XT) .ST.GENAU) GOTO 1] 
RETURN 
3 RETURN 1 
END 


### FEHLERAUSGABE - UNTERPROGRAMME ##* 


SURROUTINE ERRORL(®) 
wWRITE(6»10)9) 
100 FORMAT(IH „10X, !##=FEHLER®## WINKEL NICHT ZULAESSIG®) 
RETURN 1 
END 


SURROUTINE ERROR2(*) 

WRITE(6»100) | 
19n FORMAT(IH »10X, !##=FEHLER#F## CPI IST GROESSER ALS 16/27) 

RETURN 1 

END 


SUBRUUTINE ERRYRZ3(*) 
WRITE(6»100) 
100 FORMAT(IH »10X, !#*u=FEHLER### WIRKUNGSGRAD IST GROESSER ALS I?) 
RETURN 1 
END 


SUSROUTINE ERRURAL®) 
YRITE(&»190) 
100 FORMAT(IH »10X,. !#==FEHLER==®# CP ODER SIRKUNGSGRAD IST NEGATIV?) 
RETURN 1 
END 


SUFROUTINF ERRORS CH) 
wPITE(A»109) 
10n FNnamaT(lH „lOK.tearFEHLERere® CA IST NEGATIV 
RETURN 1 
END 


SURROUTINF ERRARSLAH®) 
CHARACTER A#13 
sRITElSs Iron) A 
19n FORMAT (///// 11m »10Xst#=n# EINGABEFEHLER #20 /1H »10XsAl8/) 


RETURM 1 
EEE EEE TEE gr ee Eee 2 222 2 202 27202202 202 202202202 
a4 ir 
##%# HINWEISE FUER PROGRAMM-SENUTZER ### 
#44 4% 


Kunsuugsanutengnenunengiihiegen et 


LAMMIN MINIMALE SCHANELLAUFZAHL AN DER FLUEGELSPITZE ; 
LAMMAX MAXIMALE SCHNELLAUFZAHL AN DER FLUEGELSPITZE 
DELLAM SCHPITTWEITE FUER DIE SCHNELLAUFZAHL AN DER FLUEGELSPITZE 


ZUR AUSGABE VON TABELLEN. IST EIN PLOTT VERLANGT (OBPL= ON?) 
"IRD DELLAM VOM PROGRAMM SO GEWAEHLT, DASS FUER JEDE DRUCK- 


a8. 


3170=C 
3180=C 
3190=C 
3200=C 
3210=C 
3220=C 
3230=C 
3240=C 
3250=C 
3260=C 
3270=C 
3280=C 
3290=C 
3300=C 
3310=C 
3320=C 
3330=C 
3340=C 
3350=C 
3360=C 
3370=C 
3380=C 
3390=C 
3400=C 
3410=C 
3420=C 
3430=C 
3440=C 
3450=C 
3460=C 
3470=C 
3480=C 
3490=C 
3500=C 
3510=C 
3520=C 
3530=C 
3540=C 
3550=C 
3560=C 
3570=C 
3530=C 
3590=C 
3600=C 
3610=C 
3A20=C 
3630=C 
3640=C 
3650=C 
3660=C 
3670=C 
I3680=C 
3690 =C 
3700=C 
3710=C 
3720=C 
3730=C 
3740=C 
3750=C 
3760=C 
37T0=C 
3780=C 
3790=C 
3800=C 
3al0=C 
3820= 


KAP 
TAU 
DELTHE 


L 

LMIN 

TT 
TLMIN 
BETT 
BETLMIN 
Z 

DELR 
GENAU 
KTR 


ORPL 


WASPL 


PARAM 


PAWERT 


PaaANZ 


CA 


EPS 


SPALTE EIN WERT BERECHNET WIRD. 

KONUSWINKEL 

ACHSWINKEL 

SCHRITTWEITE FUER DIE INTEGRATION VEBER DEM UMFANG. 

DELTHE WIRD VOM PROGRAMM SO ANGEPASST, DASS ES 360 GRAD 
OHNE REST TEILT. DIES IST DIE VORAUSSSETZUNG FUER RICHTIGE 
WIRKUNGSGRADBERECHNUNG. 


FLUFGELLAENGE: NARE - BLATTSSPITZE *] | 
LAENGE: NARBE - BLATTWURZEL »] 
FLUEGELTIEFE AN DER BLATTSPITZE =] *2 
FLUEGELTIEFE AN DER BLATTWURZEL #1 ®2 
BLATTWINKEL AN IER BLATTSPITZE #2 
BLATTWINKEL AN DER BLATTWURZEL #2 


FLUEGEL- (BLATT-) Zanl 
SCHRITTWEITE FUER DIE INTEGRATION VEBER DEM RADIUS 
GENAUISKEIT DER ITERATIVEN BERECHNUNG VON XI 
KTR=tOFfFe: 
- VERGRUESSERTER ANSTROEMWINKEL IST ZULAESSIG 
(2#ALPMIN-90 .LE.ALP-.AND.ALP.LE.180-ALPMIN) 
- ETAP=] WIRD FUER ETAP.GT.1 GESETZT 
KTR=t0oN +: 
= NUR NORMALER ANSTROEMWINKEL IST ZULAESSIG 
(ALPMIN.LEsALP.AND.ALP.LE.90) 
- FTaP.65T.ı IST FEHLER 
OBPL=MON *: DAS DIAGRAMM MIT WIRKUNGGRADKURVEN „IRD 
AUSGEGEBEN. 
OBPL='OFFt!: DAS DIAGRAMM MIT WIRKUNGSGRADKURVEN WIRD 
NICHT AUGEGEBEN. 
WASPL BESTEHT AUS 4 STELLEN: ‘ 
le STELLE=]t:! CPGES WIRD AUSGEGEBEN 
2. STELLE=:}1F: CPI WIRD AUSGEGEBEN 
3. STELLE=']Jt: ETAP KIRD AUSGEGEBEN 
4. STELLE=st1t: ETAZ WIRD AUSGEGEBEN 
REI !0': KEINE AUSGABE DES BETREFFENDEN WIRKUNGSGRADES. 
ES SIND SELIEBIGE KOMBINATIONEN "!L*tst0ü®r MOEGLICH. 
ES BESTE-T DIE MOESLICHKEIT 9 VERSCHIEDENE PARAMETER 
ANZUWAEHLEN. IN EINEM JOB WIRD DANN FUER DEN GEWAEHLTEN 
PARAMETER EINE TA3ELLE FLER JEDEN DER ANGEVEBENEN WERTE 
BERECHNET. ANSCHLIESSEND WERDEN» FALLS ANGEWAEHLT» DIE 
WIRKUNGSSRADE FUE>2 JEDEN DIESER AERTE IN E IN DIAGRAMM 
ZUM VERGLEICH SEPLOTTET. 
PARA“z= OFF 2: YIE BERECHNUNG AalIRD MIT DEN WERTEN AUS 
zEILE 200 3I 279 DURCHGEFUEHRT. 
I NER VERAENDERTE PARAMETER IST Kap 
DER VERAENDERTE PARAMETER IST L 
ARAS=ILYTN SER VERAENDERTE PARAMETER IST LMIN 
PARA“=TT DER VERAENDERTE PARAMETER IST TT 


PARAM=tKA? ' 
t 
! 
’ . 
PARA=TILMTY t2 SER VERAENDERTE PARAMETER IST TLMIN 
Gi 
t 
’ 


PARAN='L 


PARAH=IBETT NER VERAENDERTE PARAMETER IST BETT 
PARa“=ıBFTLAIN DER VERAENDERTE PAraMETER IST BETLMIN 
PARA"=ıZ ! DER VERAENDERTE ParRamETER IST Z 

PAWEPT IST EIN EINDIMENSIONALES ZAHLENFELD DER SROESSE 100. 
ES MUESSEN IMMER 190 WERTE FUER DEN PARAMETER ANGEGEBEN 
WERNEN. SEISPIEL FUER PARAM=KOFF eG 

290 DATA PAXERT /100#097 

PaANZ GI3T DIE ANZAHL DER SIGNIFIKANTEN WERTE VON PAWERT AN. 
WURDE KEIN PARAMETER GEWAEHLT, WIRD PAANZ VOM PROGRAMM 

GLEICH 1 GESETZT. 

CA IST EIN FEINDIMENSIONALES ZAHLENFELD FÜER CAC(ALP) MIT DER 
GROFESSE 112. CA : AUFTRIEBSBEIWERT. CA(ALP) ENTNIMMT MAN 
PROFILHANIDBUECHERN. FUER JEUVEN GANZEN WINKEL ALP VON =-20 BIS! 
91 GRAD „IRD DER ZUGEHOERIGE WERT FUER CA IN ZEILE 310/320 | 
ABGESPEICHERT. (SIEHE ALPMIN) I. 
EPS IST FIN ZAHLENFELD» STRUKTURIERT WIE CA» JEDICH FÜUER DIE. 
PROFILGLEITZAHL EPs BESTIMMT, 


39 | 


3830=C ALPMIN ALPMIN IST DER WINKEL» BEI DEM CA NEGATIV wIRD. 


3840=C ABGESPEICHERTE WERTE FVER CA UND EPS MIT ALP.LT.ALPMIN | 
3850=C MUESSEN VORHANDEN SEIN» IHRE GROESSE IST ABER NICHT RELEVANT 
3860=C LAMT SCHNELLAUFZAHL AN DER FLULEGELSPITZE 

3870=C CPI IDEALER LEISTUNGSSEIWERT 

3880=C ETAP PROFILREIBUNGS-VERLUSTFAKTOR 

3890=C ETA PRODUKT AUS CPI UND ETAP 

3900=C ETAZ BLATTZAHLVERLUSTFAKTOR 

3910=C QsUsTsRMIN SIEHE PROGRAMM: ### VORBEREITUNG “or 

3920=C MARKEI PROGRAMMKONTROLLMARKE. MARKEIl=1 BEDEUTET: 

3930=C VERGROESSERTER ANSTROEMBEREICH WURDE IN ANSPRUCH GENOMMEN. 
3940=C (SIEHE «TR) 

3950=C MARKE2 PROGRAMMKONTROLLMARKE. MARKE2=2 BEDEUTET: 

3960=C BEI ETAP.GT.1 WURDE ETAP=I GESETZT. (SIEHE KTR) 

3970=C MARKES3 PROGRAMMKONTROLLMARKE. MAARKE3=3 BEDEUTET: 

3980=C CAM IST ETwA 0 UND DAHER XI=1 GESETZT WORDEN. 

3990=C MARKE4 PROGRAMMKONTROLLMARKE» MARKE4=4 BEDEUTET: 

4000=C XI.LT.0 UND DAHER XI=0 GESETZT WORDEN. 

4010=C THE UMFANGSWINKEL 

4020=C THEM UMFANGSSAINKEL IN DER MITTE DES JEWEILIGEN FLUEGELELEMENTES 
4030=C RM RADINIS IN DER MITTE DES JEWEILIGEN FLUEGELELEMENTES 

4040=C LAMO EFFEKTIVE SCHNELLAUFZAHL 

4050=C AL? WINDANGRIFFSWINKEL 

4060=C SI5 DRALLFAKTOR 

4070=C DELAK FLAECHE DES FLUEGELELEMENTES 

4080=C Ak GESAMTFLAECHE BES MANTELS DES KEGELSTUMPFES 

4090=C 

4100=C 

4+110=C ANMERKUNGEN: 

4120=C 


4130=C MIT DEM PROSRAMM SuLL INTERAKTIV GEARBEITET WERDEN. DIE EINGABE DER 
4140=C GEOMETRIE UND DER RERECHNUNGSANWEISUNGEN ERFOLGT IM EDITOR: 


4150=C 1) LIST»200,300 

4160=C 2) In DEN DATA-ZEILEN NEUE WERTE EINGEBEN 

4#170=C 3) SEQUENTIELLE DATEI ERSTELLEN 

4180=C 4) geoonnnne BATCH»LDN» INPUT, AINE Pur we 

4190=C 

4200=C #1: DIE OIMENSTONSSEFAFTETEN SROESSEN L» LeMIns TT» TLMIN MUESSEN SO 
4210=C EINGESFBEN WERDEN, DASS IHRE VERHAELTNISSE DEM WINDRAD ENT- 
4220=C SPRECHEN. NIE A3SYOLUTEN BROESSEN SIND OHNE BEDEUTUNG. 
4230=C #2: FUER DEN FLUESZLTIEFENVERLAUF UND DIE SCHRAENKUNG WIRD EINE 
424U=C LINEARE VERAENDERUNG VESER DEM RAULUS ANGENOMMEN. 

4250=C 

4260=C 


42T0= END 


zn 


10.3 Eingabewerte für die Beispielplotts zum Seminarvortr und 
Tortstaberten 


E35 wurde mit dem symmetrischen Profil NACA 0015 gerechnet. Die Eingabe- 
werte für die Felder CA und EPS sind den Programmzeilen 310 bis 340 zu 
entnehmen. /4/. 


KAP TAU L LMIN TT TLMIN BETT BETLMIN Z 


45.000 45.000 10.000 3.000 .500 .500 0.000 0.000 2.000 


LAMT cPI ETAP ETA ETAZ CPGES 


1.000 .181 .137 .063 .259 .016 
2.000 .eöR .195 .085 .475 .041 
3.000 .329 293 .122 .609 .074 
4.000 .392 341 .154 .691 «106 
5.000 ..438 .380 .1838 .746 .140 
6.000 .463 .364 .191 .785 .150 
7.000 .477 .350 .192 .813 .156 
8.000 .485 313 .176 .835 «147 
9.000 .4B5 .280 .159 .852 «136 
10.000 .483 .237 .135 .866 .117 
11.000 .478 1917 .107 .878 .094 
12.000 .472 „164 .039 .888 .079 
13.000 .465 „124 .056 .896 .059 
14+.000 «460 .097 .050 .903 .045 
15.000 .455 .n31 .046 .910 .042 
16.000 .455 Re .036 .915 .033 
17.000 .4+55 .075 .03& .920 .035 
13.000 .455 .067 «033 .9.4 .031 
19.0009 .455 065 .932 «428 .030 
209.000 „455 .070 9235 «952 .033 


.050090 


DELR : 

GENAU : .0100009 
DELTHE : 22.000000 
KTR : OFF 
PARAMETER: OFF 

WERTE = 0.000 


KAP TAU L 
90.000 0.0009 10.000 
LAMT CPI ETAP 
1.000 .015 .002 
2.000 041 „014 
3.000.081 .191 
4.000  .163 „363 
5.000 .255 „528 
6.000 .327 „632 
7.000 .384 686 
8.000 .430 „709 
9.000 466 692 
10.000 493 664 
11.000 .511 .629 
12.000 .524 „574 
13.000 .533 „512 
14.000 .537 440 
15.000 .533 „354 
16.000 .534  .276 
17.000 .528  .213 
18.000 : .5138  .164 
19.000 .505 .123 
20.000 .490  .092 
DELR : .050900 
GENAU : .010000 
DELTHE  : 360.000000 

KTR : OFF 
PARAMETER: OFF 
WERTE 0.009 


LMIN 


3.000 


ETA 


.000 
.001 
.021 
.077 
.166 
.240 
.297 
.336 
.355 
.360 
.353 
.331 
.300 
.259 
.209 
.162 
«123 
.093 
.068 
.049 


TT 


.500 


ETAZ 


.259 
.475 
.609 
.691 
.746 
.785 
.813 
.835 
.852 
.866 
.878 
.888 
.896 
«903 
«910 
.915 
.920 
.924 
.928 
.932 


TLMIN 


.500 


CPGES 


.000 
.000 
.013 
.053 
.124 
.188 
.241 
.cel 
.302 
.312 
310 
.294 
.268 
.234 
.190 
.148 
„113 
.086 
.063 
.046 


BETT 


0.000 


BETLMIN 


0.000 


Z 


2.000 


52 


KAP 


LAMT 


1.000 
2.000 
3.000 
4.000 
5.000 
6.000 
7.000 
8.000 
9.000 
10.000 
11.000 
12.000 
13.000 
14.000 
15.000 
16.000 
17.000 
18.000 
19.000 
20.000 


DELR 
GENAU 
DELTHE 


KTR 


PARAMETER: 


TAU L LMIN 
50.000 60.009 10.000 3.000 
EcPI ETAP ETA 
.133 .159 .067 
.cä0 .203 .086 
.333 .306 .120 
394 .354 .155 
.437 S3t4 .179 
.462 .349 «180 
.476 315 .1l71 
.485 .2809 .157 
.487 .243 .139 
.485 .197 .113 
.48B2 1892 .0837 
.476 .114 .064 
.471 .088 .048 
.465 .066 .035 
460 . 054 .0z8 
.456 .045 .022 
.455 .034 .017 
.455 .031 .016 
.455 .034 017 
: .050000 
.: «010000 
:..20.9090°%9 
: OFF 
Tat! 
: 69.0009 30 


WERTE 


TT 


.500 


ETAZ 


.259 
.475 
.609 
.691 
.746 
.785 
.813 
«835 
.852 
.866 
.878 
.888 
.896 
.903 
«910 
.915 
«920 
.924 
.928 
.932 


.000 


TLMIN 


.500 


CPGES 


.017 
.041 
.073 
.107 
.134 
«el4l 
.139 
.1l3l 
.118 
.098 
.077 
.057 
.043 
.032 
.025 
.021 
.016 
.014 
.014 
«016 


0.000 


BETT 


0.000 


BETLMIN 


0.000 


Z 


2.000 


33 


KAP 


TAU 


L 


60.n00 30.000 10.000 


LAMT 


1.000 
2.000 
3.000 
4.000 
5.000 
6.000 
7.000 
8.000 
9.000 
10.000 
11.000 
12.000 
13.000 
14.000 
15.000 
16.000 
17.000 
18.000 
19.000 
20.000 


CPI 


.032 
.ce]l 
.319 
.4+00 
.454 
.483 
.506 
.511 
+.508 
.504 
.498 
.492 
.«4R5 
.477 
.4+70 
.463 
.458 
.455 
.455 


ETAP 


.007 


LMIN 


3.000 


ETA 


.0Cc0 
.036 
.l27 
.222 
.292 
.324 
.295 
.259 
.221 
.186 
.154 
.lel 
.091 
.068 
.048 
038 
.027 
.0285 
.031 


TT 


.500 


ETAZ 


«259 
.475 
«609 
.691 
+746 
.785 
«813 
«835 
.852 
.866 
.878 
.888 
«896 
.903 
-.910 
.915 
.920 
.924 
.928 
.932 


TLMIN 


.500 


CPGES 


.000 
.017 
.077 
.154 
.218 
.254 
.261 
.246 
.220 
.192 
.163 
.137 
.109 
.082 
.062 
.044 
.035 
.025 
.026 
.029 


BETT 


0.000 


BETLMIN 


0.000 


zZ 


2.000 


34 


KAP 


60.000 


LAMT 


1.000 
2.000 
3.000 
4.000 
5.000 
6.000 
7.000 
8.000 
9.000 
10.000 
11.000 
12.000 
13.000 
14.000 
15.000 
16.000 
17.000 
18.000 
192.000 
20.000 


TAU 


0.000 


.022 
.060 
.1283 
.253 
.343 
.417 
.468 
.503 
.524 
534 
.537 
.534 
.524 
.509 
.489 
.4+65 
.454 
+.454 
.454 
.454 


lE 


10.000 


LMIN 


3.000 


ETA 


.000 
.00l 
.050 
.152 
«250 
«323 
. 364 
.377 
.366 
«336 
.289 
.229 
.164 
“114 
.077 
.049 
.040 
«043 
«030 
.020 


Tr 


.500 


ETAZ 


.259 
.475 
«609 
.6891 
.746 
.785 
.813 
.835 
.852 
.866 
.878 
.888 
.896 
.903 
«910 
.915 
.920 
.924 
.928 
.932 


TLMIN 


.500 


CPGES 


.000 
.000 
.031 
.105 
.186 
.253 
.296 
3:15 
.3l2 
.291 
.253 
.203 
.147 
.103 
.070 
.045 
.037 
.040 
.028 
.019 


BETT 


0.000 


BETLMIN 


0.000 


zZ 


2.000 


35 


KAP 


TAU LE 


90.000 45.000 10.000 


LAMT CPI 
1.000 +.030 
2.000 .086 
3.000 .192 
4.000 .301 
5.000 .332 
6.000 .441 
7.000 .482 
8.000 .511 
9.000 .523 
10.000 .535 
11.000 .53& 
12.000 .526 
13.000 .511 
1+.000 .492 
15.000 .471 
16.000 .4öR 
17.000 .455 
13.000 .455 
19.000 .455 
20.000 .455 
DELR : 
GENAU : 
DELTHE 29 
KTR : OFF 
PARAMETER: Kar 
WERTE : 


ETAP 


«009 
.184 
.4+19 
.532 
.666 
.683 
.659 
.603 
.522 
.420 
.310 
.2l7 
.150 
.101 
.071 
.053 
.038 
.040 
.039 
«032 


«050000 
.N10000 
.0000099 


90,299 


LMIN TT 
3.000 «500 
ETA ETAZ 

.000 .259 
«030 .475 
.114 .609 
.216 .691 
.294 .746 
.337 .785 
.351 .813 
.338 .835 
.302 .852 
.246 .866 
151 .878 
.126 .838 
.035 «896 
. 055 «903 
.037 .910 
.027 .915 
.019 .920 
.020 .924 
.019 .928 
.016 «932 
67.500 


TLMIN 


.500 


CPGES 


.000 
.014 
.069 
.149 
.219 
.264 
.285 
.282 
.257 
.213 
.159 
.lı2 
.076 
.050 
.034 
.025 
.017 
.019 
.013 
«015 


45.000 


BETT 


0.000 


BETLMIN 


0.000 


zZ 


2.000 


36 


KAP 


AT.Sou 


LAMT 


1.000 
2.000 
3.000 
4.000 
5.000 
h.Nn00 
7.000 
3.000 
9.900 
10.900 
11.000 
12.000 
13.000 
14.200 
15.000 
15.200 
17.000 
12.000 
20.000 


Ta’ 


+5.090° 


L 


19.090 


ETAP 


.072 
267 
.463 
Sie 
.5399 
oe 
>03 
.432 
Bla 
EN, 
.260 
205 
.11# 
.985 
.075 
.063 
.059 
.052 
.145 


LAT 


3.120 


ETA 


.)99 
«Dal 
«1.72 
.2x7 
.2e6 
«256 
.239 
.c2ıl 
«1-2? 
.150 
el] 
«0-2 
.Nn4 
.0a7 
lan 
«032 
.030 
«02H 
.072 


.c59 
.4/5 
.5093 
.7465 
.755 
.813 
.535 
.öb2 
.300 
.873 
«910 
.915 
.9c0 
.«9L# 
.3c9 
932 


TLHIN 


«>99 


CPGES 


.002 
«0 Ic 
«105 
171 
.el4 
.cL+ 
.21n 
„1928 
13] 
079 
«058 
«042 
.037 
030 
.023 
.024 
.Dcl 


gETT 


0.090 


SETLMIM 


0.000 


Z 


2.00 U 


Sr 


= Ir 2 
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